radiation->correct();
radiationF = radiation->solidSh();

